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I.  INTRODUCTION 


We  consider  two  nearly  identical  Fabry-Perot  laser  oscillators,  which 
are  coupled  by  partial  injection  of  the  output  of  one  laser  into  the  second 
and  vice  versa.  Spencer  and  Lamb^  have  shown  that,  under  some  circumstances, 
a  steady-state  solution  can  be  obtained  where  the  two  lasers  are  locked  to  a 
fixed  frequency  and  phase  relation.  This  technique  can  then  be  used  to  obtain 
the  coherent  combination  of  the  output  from  the  two  lasers. 

In  the  theory  of  Spencer  and  Lamb,  it  is  assumed  that  there  is  a  fixed 
phase  difference  (herein  denoted  ip  =  -  tt / 2 )  between  the  transmitted  and 
reflected  radiation  at  the  coupling  mirrors.  This  boundary  condition  is 
appropriate  when  the  coupling  mirrors  are  equivalent  to  a  lossless  Fabry-Perot 
etalon  (see  Appendix).  A  more  general  treatment  is  to  consider  the  value 
of  ij>  to  be  arbitrary.  This  approach  is  used  herein. 

The  equations  that  describe  the  performance  of  coupled  laser  oscillators 
operating  on  a  single  frequency  are  first  noted.  Weak  coupling  is  assumed.  A 
third-order  saturation  model2  is  used  for  both  homogeneous  and  inhomogeneous 
laser  media.  Analytic  and  numerical  solutions  are  then  obtained  for  special 
cases  with  emphasis  on  the  stability  range  for  steady-state  operation. 

Finally,  the  stability  of  multiline  coupled  oscillators  is  discussed. 


[I.  THEORY 


Equations  that  define  the  performance  of  single  line  coupled  resonators 
are  deduced  herein. 


A.  RESONATOR  PROPERTIES 


Coupled  resonators  are  illustrated  in  Fig.  la.  The  coupling  occurs 
through  mirrors  and  M21’  The  e^-ectr^c  field  at  these  mirrors  is  noted  in 

iij>  j 

Fig.  lb.  The  reflected  electric  field  at  and  ^21  *s  denoted  E  e  f°r 

1 4>  2 

and  E^e  for  ^l*  The  corresponding  transmitted  electric  field  in  Mj2 
1  /2  i  (4> )  i  u  i  ( «t> .  +U> ) 

and  ^  is  (T/R)  E^e  and  (T/R)1  E^e  .  Here,  T  and  R  repre¬ 
sent  net  intensity  transmission  and  reflection  coefficients  while  ^  — d, r 

is  the  net  phase  difference  between  the  transmitted  and  reflected  waves. 

Values  of  T,  R,  and  t  are  given  in  the  Appendix  for  cases  where  mirrors  Mj2 
and  M21  form  a  Fabry-Perot  etalon.  When  there  are  no  losses,  T  +  R  =  1 
and  t  =  ±  tt/2  .  The  latter  result  is  consistent  with  assumptions  in  Ref.  I. 
However,  when  losses  occur,  T  +  R  <  1  and  t  depends  on  the  optical  path  length 
2irh/X  (Appendix).  Hence,  consideration  is  given  herein  of  cases  for  which 
|  p  |  4  tt/2  . 

Vector  addition  of  the  reflected  and  transmitted  wave  in  each  resonator 


(Fig.  2)  provides  the  net  perturbation,  per  round  trip,  caused  by  the 
injection  process.  For  weak  coupling,  (T/R)*^  <<  1,  the  result  is 


1/2 


AEj  =  (T/R)  E2  cos(4>  +  f ) 


(la) 


1  /2 

At,  =  (T/R)  (E  /E  )  sin(t  +  t ) 


(lb) 


% 


RESONATOR  1  RESONATOR  2 


AE2  =  (T/K)1/2E1  cos  (<{,  -  ip) 


1  /2 

A<p 2  =  -(T/R)  (e1/e2)  sln  ~ 


(lc) 


(Id) 


where  <j>  =  $2  ~  <t>^«  The  round-trip  time  in  each  resonator  is 


tr.n  =  2  Ln/c 


where  n  =  i,2.  The  empty  cavity  resonance  freqency  for  each  resonator  is 


(2) 


=  ( 2ttN l  -  *  )(c/2L  ) 


n2  =  (2ttN2  -  $2)(c/2L2) 


(3a) 


(3b) 


where  and  N2  are  integers.  Here,  $ ^  and  $2  denote  the  net  phase  change 
induced  by  reflections  during  a  round  trip  within  resonators  1  and  2, 
respectively.  Both  homogeneous  and  inhomogeneous  media  are  considered.  For 
Later  reference,  we  ..itroduce  the  normalized  frequency  relations 


K  =  (w-v)/y 


A  =  (u)-Q  )/y 
n  n 


*  A2  -  Aj  =  (Q1  -  n2)/y 


(4a) 


(4b) 


(4c) 


-Here  a),  v,  and  y  denote  line  center  frequency,  laser  frequency,  and  homogen¬ 
ous  Line  width  (HWHM),  respectively. 
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B. 


UNSTEADY  RESONATOR  EQUATIONS 


Single-mode  operation  of  a  Fabry-Perot  resonator  has  been  described  by 
Lamb.2  These  equations  can  be  generalized  to  the  case  of  two  coupled  resona¬ 
tors  by  inclusion  of  Eqs.  (la)  through  (Id).  The  result  is 


i 

— -  =  (a.  -  8,1.)  E  +  M  E  cosU  +  ij>) 
dt  1  11  1.1/ 

dE? 

aT  '  U2  -  e2V  L  *  M2E!  CO,(*  "  *> 


(5a) 


(5b) 


d  <J> 

=  y(C  -  Aj )  +  aj  ~  Pjlj  +  M1(E2/E1 )  sin(*  +  *) 


(5c) 


d<J>_ 

dT  “  ‘  V  +  °2 


P2i2  "  M2(R1/E2)  sin(<t> 


*> 


( 5d ) 


^  +  a2  -  +  Pjlj  -  P2I2  -  M2(E1/E2)  s i n ( -  i|i) 


M  (E  /Et )  sin(  <f»  +  #) 


(5e) 


where  I 


En2  and 


M  =  (T/R) 1 /2/t  (5) 
n  t  ,n 

The  last  term  on  the  right  hand  side  of  Eqs.  (5a)  through  (5d)  is  deduced  by 
dividing  Eqs.  (la)  through  (Id)  by  tr>n  to  obtain  the  rate  of  change  per  unit 
t  ime. 
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The  quantities  a  ,  8  a  ,  p  are  functions  of  F  and  resonator  parameters 
1  n*  n  n  n 

F,  ,  F0  ,  and  F,  as  defined  in  Table  I.  When  i|>  =  -  ir/2,  Eqs.  (5a)  through 
l ,n  *  4,n ’  i,n 

(5e)  agree  with  the  corresponding  equations  in  Ref.  1  except  for  the  satura¬ 
tion  model,  which  herein  is  limited  to  the  third-order  (in  E)  theory  of 
Lamb.  The  third-order  theory  is  accurate4  for  1  <  Fj  n/F2  n  ~  (denoted 
N/Nt  in  Ref.  2),  but  is  expected  to  give  qualitatively  correct  results  for 
larger  values  of  this  parameter. 

C.  STEADY-STATE  SOLUTION 

A  steady-state  solution  of  Eqs.  (5a)  through  (5e)  is  obtained  by 
assuming  that  the  time  derivative  terms  are  zero.  The  resulting  equations  can 
be  put  in  the  form 


1  /2 

ll  02tal  +  VW  cos +  ♦>! 

_  i/o  (2a) 

2  Bj[a2  +  M2(I1/I2)  cos(<j>  -  *)] 

1  /2 

yA.  -  M.  (I  /I  )  sin($  +  <|>) 

C - ! - L_2_! - m - (7W 

Y  +  ((Jj/C)  -  (  pj  /  C6 1  )  f a  i  +  M1(I2/I1)  cos(<J>  +  ij>)] 


1/2 


Ij  =  [aj  +  Mj(I2/I1)  cos(<f>  +  !|))1/8j 


(7c) 


l2  (I2/T1)  l\ 


(7d) 
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TABLE  I. 

Coefficients  that  define 
from  Ref.  2 

gain  medium  saturation 

Coefficient 

Homogeneous 

Medium 

Inhomogeneous 

Medium 

an 

®>Fl ,n  -  F2,„ 

-,2 

e  C  F,  -  F. 

1 ,n  2,n 

B 

^2f3  n 

(1  +&)  f3 

n 

an 

j  ,n 

«*F1,» 

•J  y  II 

2  Fl,n  D(C) 

pn 

?6n 

5^?F_ 

3,n 

(g00*/tr)n  (g  =  d£nl/dx) 

(gc^/cr 

F,  (assumes  I  normalized  by  saturation  value) 

l  ,n 

(a)  -  v)/y 
(1  +  52)"1 

d  -t2  C  x2 

Dawson  integral  =  e  ^  J  eX  dx 

o 

line  center  -  zero  power  -  gain 
threshold  gain  =  -  £n(R^R)^2/£ 
length  of  gain  region 


(7e) 


yA  =  a2  -  <jj  +  p  j  I  j  -  p9I2  "  s  i  n  ( <t>  -  i|») 

-  M1(I2/I1)1/2  sin(  (j>  +  *) 

A2  =  A  +  Aj  ( 7  f ) 

The  procedure  for  obtaining  a  solution  Is  as  follows. 

The  quantities  $,  <J/,  M,  ,  M  ,  A.,  y,  F.  ,  F.  and  F9  „  are  specified. 

1  l  L  1  ,n  z  ,n  5  >" 

For  an  inhomogeneous  medium,  the  doppler  half  width  y  is  also  specified.  The 

d 

value  of  5  is  estimated  and  T ^ / 1 2  is  obtained  from  Eq.  (7a).  A  new  estimate 
for  £  is  then  obtained  from  Eq.  (7b).  The  procedure  is  repeated  until  the 
estimate  for  £  converges.  The  quantities  Tj,  I2,  A,  and  A2  are  then  found 
from  Eqs .  (7c)  through  (7f).  The  variables  1J1,  Mj,  Aj,  Aj,  y,  Fj  n,  F2  n, 
and  F^  n  are  resonator  and  gain  medium  properties.  The  above  procedure 

provides  the  corresponding  values  of  laser  frequency  v,  phase  angle  <J> ,  and 

laser  intensities  Ij  and  I2<  Equation  (7a)  is  fourth  order  in  (Ij/I2)'^  so 
that  four  steady-state  solutions  for  (Ij/l^)1^2  are  possible,  in  principle. 
However,  in  all  cases  considered  herein,  no  more  than  one  stable  steady-state 
solution  has  been  found  for  a  given  set  of  initial  conditions.  Stability  is 
examined  in  the  next  section. 

D.  STABILITY  CRITERION 

The  stability  of  the  steady-state  solutions  can  be  determined  by 
considering  the  effect  of  small  perturbations. 

Let  E|  =  E10  +  AEj,  E2  =  E2q  +  AE2 ,  and  $  =  $  +  A$  where  subscript  zero 

denotes  the  steady-state  solution  and  A  ()  denotes  an  unsteady  perturbation. 
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In  the  limit  <|>  =  -ir/2,  these  equations  agree  with  corresponding  equations  In 
Ref.  1.  Perturbations  are  damped,  and  the  system  Is  stable,  If  the  roots  of 
the  characteristic  equation 


a  .  ~  6  r 
1  ij  ij 


(10) 


have  negative  real  parts. 

E.  SIMILITUDE 

The  steady-state  equations  depend  on  $  and  ip  through  terms  of  the  form 
cos(<|>  ±  i|i)  and  sin(if>  ±  i|>).  Equivalent  cases  are  now  noted.  Let  f ( <J> ,  <J<) 
denote  a  dependent  variable.  Equations  (7a)  through  (7f)  Indicate 

f  (  ±  it,  ij;  ±  ir)  =  f  ( 4> ,  Ii>)  (11) 


For  the  special  case  A^  =  0, 

f(-$,  “♦)  -  f(<>,  <(>)  f  =  T1  ’  l2 

=  -f (♦,  i|>)  f  -  ?,  A,  J 


(12a) 

(12b) 


where  $  is  defined  in  Eq.  (IS). 

F.  BEAM  COMBINATION 

The  output  from  each  of  the  coupled  resonators  can  be  directed  to  anten¬ 
nae  which  are  either  separated  (Fig.  3a)  or  collocated  (Fig.  3b).  When  a 
coupled  resonator  is  operating  under  steady  state  conditions  (l.e.,  lock  in), 


V  1 


the  configuration  in  Fig.  3a  constitutes  a  phased  array  while  the  configura¬ 
tion  in  Fig.  3b  constitutes  a  technique  for  coherent  beam  combination. 
Factors  affecting  the  beam  quality  of  the  latter  configuration  are  discussed 
herein  for  both  single  line  and  multiline  coupled  resonators. 

The  beam  quality  of  the  output  beam  indicated  in  Fig.  3b  can  be 
expressed 


-  ^ideal^ ^actual 


where  *actuaj_  is  far-field  center  line  intensity  and  is  the  value  for  a 

combined  beam  with  uniform  near-field  phase  and  intensity.  Values  of  beam 
quality  above  1  indicate  a  lack  of  near-field  phase  uniformity.  Assume  that 
the  output  beams  from  the  resonators  in  Fig.  3b  are  equal  in  area  and  inten¬ 
sity  and  illuminate  subappertures  (A1  and  A2)  of  equal  area.  Let  <j>A  and  $ 

1  A2 

denote  the  phase  of  the  radiation  at  Aj  and  A2,  respectively,  and  assume 

2  q 

(<{>*  -  <{>A  )  «  1.  The  beam  quality  of  the  combined  beam  is  thenJ 

2  A1 


+  °^A  "  $A  * 

A2  A1 


BQ  -  1  =  - 


Equation  ( 14)  can  be  evaluated  as  a  function  of  coupled  resonator  geometry  and 


gain  medium  properties. 


w 


For  a  single  line  stable  coupled  resonator  the  phase  of  the  output  beam 
at  stations  and  can  be  shown  to  be  <J>^  +  (Atf^/2)  and  <J> ^  +  (44^/2)  -  ki,  ^ , 

respectively,  where  k  =  2tt / A .  Similarly,  the  phase  of  the  output  beam  at  M22 
and  A2  is  <j>2  +  ( A<J> 2 /2 )  and  $2  +  ( A(f>2 /2 >  “  resPectl vely*  The  phase 

difference  between  the  radiation  at  A2  and  Aj  is  then 

*A2  "  *A  =  ♦  +  ^R^1/2  *  -  k(*2  "  £1}  °5a) 

where  <j>  =  4>  2  -  <j>  ^  and 

A$^  —  A$ .  1^1  ^7 

<t>  5  - r-pr  =  -  «•  [•=—  sin(<J>  -  ijO  +  —  sin(4>  +  i|»)]  (15b) 

2(T/R)  1  2  1 

=  -  sin  $  sin  ip  (E ^  =  E2)  (15c) 

Substitution  of  Eq.  (15a)  into  Eq.  (14)  defines  the  beam  quality  of  the 

combined  beam.  In  a  single-line  stable  resonator,  with  4  L2,  it  can  be 

shown  that  <j>  4  0,  However,  the  condition  <t>  ~  <f>.  =0  (i.e.,  good  beam 

A2  1 

quality)  can  be  achieved  by  proper  choice  of  l 2  -  in  Eq.  (15a).  In  a 
multiline  stable  coupled  resonator,  with  Lj  4  L2,  however,  the  values  of 
k,  <J> ,  and  <j>  differ  for  each  wavelength.  Hence,  it  is  not  possible  to  adjust 
H2  -  l so  as  to  obtain  good  beam  quality  for  all  wavelengths.  Therefore,  a 
necessary  condition  to  achieve  good  beam  quality,  for  the  case  of  a  multiline 
stable  coupled  resonator,  is  to  operate  under  conditions  where  =  L2  and 
l  =  l  .  In  this  case  <J> ,  <J>  and  -  A  are  small  for  each  wavelength  and 

12  A2  Aj 

good  beam  quality  is  obtained  provided  the  solution  is  stable. 
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III.  SPECIAL  CASES 


Numerical  methods  are  generally  needed  to  evaluate  coupled  resonator 

performance.  Explicit  expressions  can  be  obtained  in  special  cases.  The 

2 

cases  cos($  +<!')  =  cos(4>  -  ip)  and  <J>  <<  1  are  considered  here. 

A.  ASSUMPTIONS 

Let  L  =  (Lj  +  L2)/2  denote  the  average  length  of  resonators  1  and  2  and 


(L2  -  L^/L  =  0(X/L)  «  1 


It  follows  that 


tr>n  =  tr  [1  +  0(X/L) ] 


Mn  =  M [ 1  +  0(X/L) ] 


where  tr  =  2L/c  and  M  =  (T/R)^2/tr.  We  also  assume  that  the  gain  medium  and 

the  mirror  configuration  for  resonators  1  and  2  are  nominally  identical.  We 

therefore  replace  M_,  a  .  6  ,  a  ,  and  p  in  Eqs.  (5)  and  (7)  by  the  mean 

n  n  n  n  n 

values  M,  a,  8,  o,  and  p.  Under  these  conditions  the  dependent  variables 


1 1  >  ^2 ’  ^  ' 


are  functions  of  the  parameters 


To  facilitate  computations,  $  is  sometimes  specified  and  the  corresponding 
value  of  Aj  (or  h  )  is  determined. 

B.  CASE  cos(4>  +  ij>)  =  cos(j)  -  ijt) 

To  simplify  Eq.  (7a),  we  assume  that  $  and  ij>  satisfy  the  relation 

cos(<{>  +  4»)  =  cos(<J>  -  =  cosiJ>  cosi|i  (19 

Equation  (19)  is  equivalent  to  siniji  sinij;  =  0  and  is  satisfied  by 


<p  -  0,  it;  iJi  =  arbitrary 


(20a 


and  by 


=  0,  tt ;  *  arbitrary 


(20b 


Equation  (20a)  includes  the  case  where  Lj  =  and  h  is  variable. 

Equation  (20b)  includes  the  case  where  L|  +  l^  and  kh  has  a  value  equal  to  an 
integral  number  of  half  wave  lengths.  Substitution  of  Eq.  (19)  into  Eq.  (7) 


yields 


where 


D  =  (M/a)  eosf  (21 

and  the  root  (Ij/I0)  '  =  -1  has  been  omitted.  Equation  (21b)  provides  real 

positive  solutions  when  -  ^  <  D  <  0.  Several  numerical  examples,  however, 
indicate  that  these  solutions  are  unstable.  We  therefore  pursue  the  case  Ij 
I9  =  I.  Substitution  of  the  latter  relation  into  Eq.  (7)  yields 

yAj  -  M  si  n(  j,  +  J/) 

’  Y  +  (a/f;)  -  [p/(CS)][a  +  M  cos(<fc  +  ij>)] 


yAj  -  M  sin($  +  lb ) 

Y  +  “  M  cos($  +  iji) 


homogeneous 
med i um 


(22 


I  =  [n  +  M  cos  ( <j>  +  ^)  ]  /  B 


(22 


YA  =  -2M  sinij)  cos>i) 


(22 


£  =  n 


(22 


Equations  (22a)  through  (22e)  provide  a  solution  for  £, 
M,  $ ,  i|>,  F  ,  and  F,  are  specified. 


and  A  when  y. 


2  7 


When  tp  =  0,  ir  and  $  is  arbitrary  the  roots  of  Kq .  (10)  are 


r^  =  — 2M  eos<|>  cos-j 


r 2  rl  "a 


r3  =  2r,  -  2a 


The  coupled  resonator  solution  is  stable  when  r.  <  0  or 


(  2  la  ) 


(lib) 


me) 


-ir/2  <  4>  <  tt/2  ip  =  0 


tt/2  <  4>  <  3tt/2  ip  =  it 


(24a) 


(24b) 


Note  the  identity 


Q2  Ql 
2 tt  (c/2L) 


L1  l2 


~  yA(T/R) 
2m  M 


It  follows,  from  Eqs.  (22d),  (24),  and  (25)  that  the  tuning  range  for  stable 
operation  is 


m  m  mi/2 

\T?  <  7T  m 


1/2 

tuning  range  increases  linearly  with  (T/R)  '  . 


Similarly,  when  <j>  =  0 ,  it  and  tp  is  arbitrary  the  roots  of  Eq .  (10)  are 


r  =  -2  (  2a  +  b)  /  3 


(27a) 


*.  \  \  •.  *.  *.  %  \  V*. 
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V- 


.  ,.2  .  .1/2 

r2  ^  =  -b  ±  (b  -  4c) 


(27b) 


where 


b  -  a  +  3M  cos<f>  cosi)/ 


(27c) 


dj 


2  2 

c  =  M  (1  +  cos  iji)  +  M  cosij)(a  cosi))  +  pi  sinijj)  (27d) 


Consideration  of  the  root  ^  indicates  that  the  system  is  stable  when  the 
conditions  b  >  0  and  c  >  0  are  satisfied.  The  condition  b  >  0  is  satisfied  by 


cosij)  >  - 


COS!))  < 


ij>  =  0 


4>  =  TT 


O 

In  the  limit  (a/3M)  «  1,  the  latter  stability  region  becomes 


1*1  <  f  *  +  “(Iff)3 


11  >  1 4)1  >  2  3M  +  0  ^  3M^ 


The  condition  c  >  0  is  satisfied  by 


( 23a) 


(28b) 


(29a) 


(29b) 


M  2 

COSi)>  COSlJ)  > - (1  +  COS  !))) 

3 


-jp-  sinij)  cosiji 
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For  a  homogeneous  medium  with  <p  =  A 


1 


0,  Eq.  (30)  becomes 


COSlp  > 


1  -  (y  +  F ^ ) C 1  +  cos  \p)/a 


Y  +  F, 


(31a 


1  -  2M  / [ a( y  +  F2)l 


Y  +  F, 


[1  +  0(- 


y  +  F, 


*  o( 


a( y  +  F2) 


n 


(31b 


When  a/3M  is  small,  values  of  tp  for  stable  operation  are  determined  by  the 

condition  b  >  0.  When  a/3M  is  large,  values  of  ip  for  stable  operation  are 

determined  by  the  condition  c  >  0.  The  range  of  ip  for  stable  operation  first 

increases  and  then  decreases  with  increase  in  a/M  (Fig.  4).  In  cases  where 
2 

5  <<  1,  the  substitution  a  =  Fj  -  F2  simplifies  the  evaluation  of  stability. 

C.  CASE  <t>2  «  1 
2 

We  assume  (p  «  1  and  express  dependent  variables  in  the  form 


E  /E  =  1  +  *E  (32a 

a  =  a  (1  +  *a  )  (32b 

8  =  S0(l  +  ♦Bjj)  (32c 


where  subscript  zero  denotes  values  for  $  =  0  and  where  a  j |  is  not  to  be 

confused  with  Eq.  (9a).  Substitution  into  Eqs.  (7a)  through  (7f)  and  neglect 

2 

of  terms  of  order  $  yields 


F. 


1  1 


-A 

1  +  2A 


t  an 


<P 


(3  3a 
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RESONATOR  SEPARATION,  (180/tt)  i,  degrees 


SATURATION  PARAMETER,  (F1  -  F2)/M 


Fig.  4. 


Effect  of  saturation  on  range  of  ip  for  stable  operation.  A 
homogeneous  medium  is  assumed  with  input  conditions  given  in  Tabl 


II. 
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A— *. 


-^7  =  E  -7-^  ( 1+A)  +  M  slrnjjl  -  M  cos\|> 
la  1  1  S  1 


<J>  =  -  <j>(cos<j>  -  E  sinij;) 


where 


A  =  (M/a^cosiJ) 


Equation  (33a)  indicates  the  departure  from  the  condition  I ^  / X 2  =  1  as  <J>  is 


increased  from  zero.  Equation  (33b)  indicates  the  relation  between  $  and  the 


tuning  parameter  A.  Additional  expressions  for  dependent  variables  can  be 


readily  obtained. 


IV.  RESULTS  AND  DISCUSSION 


The  effect  of  variations  in  $  and  \jj  on  steady-state  coupled  resonator 
performance  has  been  evaluated  for  initial  conditions  given  in  Table  TI. 
Calculations  were  conducted  at  Intervals  of  A<f>  =  Aiji  =  2°.  The  results  are 
given  in  Tables  III  and  IV  for  homogeneous  and  inhomogenous  media,  respec¬ 
tively,  and  are  discussed  herein.  The  end  of  each  stable  operating  regime, 
denoted  by  an  asterisk  in  Tables  III  and  IV,  is  correct  to  within  2°  due  to 
the  computational  interval. 

A.  <fr  FIXED  AND  VARIED 

The  present  computations  correspond  to  keeping  =  0)  and  L2*-A2~  ^ 

fixed  while  varying  h(tj/). 

The  net  output  is  proportional  to  Ij  +  I2  and  is  a  maximum  when 
4>  =  i|>  =  0.  When  <j>  =  0,  an  increase  in  |i{j|  results  in  a  decrease  in  net  output 
and  stable  operation  is  limited  to  the  range  -96°  <  <  96°.  With  Increase 

in  <J>  the  peak  output  continues  to  occur  at  ip  =  0  (except  when  <t>  =  90°)  but  is 
reduced  from  the  output  at  $  =  ijj  =  0.  The  range  of  i|j  for  stable  operation  is 
also  reduced  with  increase  in  4>. 

The  data  in  Tables  III  and  IV  correspond  to  a  fixed  degree  of  satur¬ 
ation,  F ^ / F 2 *  The  effect  of  this  parameter  on  the  range  of  \ p  for  stable 
operation  is  indicated  in  Fig.  A  where  stability  limits  are  plotted  versus  ( F ( 
-  F2)/M.  When  <j>  =  0,  the  stable  operating  regime  limit  increases  from  ij;  =  ± 
90°  at  (Fj  -  F2>/M  =  0  (threshold)  to  =  ±  130°  at  (Fj  -  F2)/M  =  2.  These 
results  can  be  obtained  by  taking  b  =  0  in  Eq .  (27c).  The  range  of  ip  for 
stable  operation  decreases  with  further  increase  in  (Fj  -  F2)/M  and  is  deduced 
by  taking  c  =  0  in  Eq.  (27d).  An  increase  in  A  results  in  a  decrease  in  the 


TABLK  It.  Input  for  numerical  computations 


= 

0.50 

X 

108 

rad/s 

Yd 

= 

1.00 

X 

109 

rad/s  (Table  IV) 

M 

= 

0.10 

X 

108 

s'1 

F1 

= 

0.20 

X 

108 

s'1 

= 

0.17 

X 

108 

-  0.60  x  108  s"1  (Fig.  4) 

F2 

= 

0.16 

X 

108 

s_1 

A. 

_ 

0 

1 
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III.  Coupled  resonator  performance  for  case  of  homogeneous 
medium.  Input  conditions  given  in  Table  II.  Asterisk 
denotes  end  of  stable  operating  regime  to  within  2°. 


(a)  4>  Fixed 

and  Varied 

180  , 
—  * 

K 

*1 

h 

1 1  +  I2 

A2 

0 

0.000 

0.700 

0 . 700 

1.400 

0.000 

±30 

TO. 087 

0.635 

0.635 

1  .270 

0.000 

±60 

TO. 142 

0.448 

0.448 

0.896 

0.000 

±90 

TO. 152 

0.186 

0.186 

0.372 

0.000 

±96* 

TO. 148 

0.132 

0.132 

0.264 

0.000 

-114* 

0.153 

0.240 

0.239 

0.479 

0.081 

-90 

0.091 

0.362 

0.156 

0.518 

-0.143 

-60 

0.053 

0.479 

0.198 

0.677 

-0.240 

-30 

0.000 

0.637 

0.486 

1.123 

-0.198 

0 

-0.087 

0.635 

0.635 

1.270 

-0 . 200 

30 

-0.165 

0.486 

0.642 

1.128 

-0.189 

60 

-0.256 

0.157 

0.448 

0.605 

-0.252 

62* 

-0.298 

0.099 

0.393 

0.492 

-0 .315 

-90* 

0.066 

0.539 

0.331 

0.870 

-0.044 

-60 

0.000 

0.497 

0.176 

0.673 

-0.291 

-30 

-0.053 

0.479 

0.198 

0.677 

-0.381 

0 

-0.142 

0.448 

0.448 

0.896 

-0.346 

30 

-0.256 

0.157 

0.448 

0.605 

-0.371 

32* 

-0.309 

0.093 

0.393 

0.486 

-0.426 

-6*/ 186* 

TO.  109 

0.231 

0.118 

0.349 

TO. 425 

-4/184 

TO. 124 

0.220 

0.146 

0.366 

TO. 41 1 

-2/182 

TO. 138 

0.205 

0.168 

0.373 

TO.  404 

0/180 

TO. 152 

0.186 

0.186 

0.37  2 

TO. 400 

2 / 1 78  * 

TO . 1 67 

0.162 

0.199 

0.361 

TO. 400 

4/176 

TO. 188 

0.131 

0.208 

0.339 

TO. 404 

TABLE  III.  Coupled  resonator  performance  for  case  of  homogeneous 
medium.  Input  conditions  given  in  Table  II.  Asterisk 
denotes  end  of  stable  operating  regime  to  within  2° 

( Continued ) . 


(b)  Fixed  and  $  Varied 


180 

r. 

[1 

l2 

Vl2 

A2 

0 

0.000 

0.700 

0.700 

1.400 

0.000 

±30 

40.087 

0.635 

0.635 

1.270 

+0.200 

±60 

40.142 

0.448 

0.448 

0.896 

+0.346 

±90* 

40.152 

0.186 

0.186 

0.372 

+0.400 

-92* 

0.096 

0.357 

0.166 

0.523 

0.376 

-90 

0.091 

0 . 3o2 

0.156 

0.518 

0.385 

-60 

0.053 

0.479 

0.198 

0.677 

0.381 

-30 

0 . 000 

0.637 

0.486 

1 .123 

0.198 

0 

-0.087 

0.635 

0.635 

1.270 

0.000 

30 

-0.165 

0.486 

0.642 

1.128 

-0.189 

60 

-0.256 

0.157 

0.448 

0.605 

-0.371 

62* 

-0.298 

0.099 

0.393 

0.492 

-0.422 

-144* 

0.164 

0.242 

0.277 

0.519 

0.134 

-140 

0.159 

0.279 

0.298 

0.577 

0.136 

-130* 

0.145 

0.360 

0.338 

0.698 

0.148 

-74* 

0.028 

0.538 

0.261 

0.799 

0.243 

-60 

0.000 

0.497 

0.176 

0.673 

0.291 

-30 

-0.053 

0.479 

0. 198 

0.677 

0.240 

0 

-0.142 

0.448 

0.448 

0.896 

0.000 

30 

-0.256 

0.157 

0.448 

0.605 

-0.252 

12 

-0.309 

0.09  3 

0.393 

0.486 

-0.335 

★  .  .  ,  * 
-b i)  /  2*40 

40 .066 

0.5  39 

0.331 

0.870 

±0.044 

-30/230 

40.081 

0.491 

0.28  3 

0.774 

±0.065 

-40/220 

40.090 

0.431 

0.224 

0.655 

±0.095 

-30/210 

40.09 l 

0.362 

0.156 

0.518 

±0 .143 

1 

r-j 

ro 

40.083 

0.304 

0.099 

0.403 

±0.213 

-1O/190 

40 . 142 

0.271 

0.232 

0.50  3 

±0.028 

o/l  80 

40. 1  52 

0.186 

0.186 

0.372 

0.000 

*  _  w 

■»  /  1  1  h 

h).  188 

0.131 

0.208 

0. 339 

+0.088 

TABLE  IV.  Coupled  resonator  performance  for  case  of  inhomogeneous 
medium.  Input  conditions  given  in  Table  II.  Asterisk 
denotes  end  of  stable  operating  regime  to  within  2°. 


(a)  4>  Fixed 

and 

Varied 

_ 

_ 

180  , 
—  * 

180  , 

__  * 

S 

*1 

l2 

M  +  I2 

A2 

0 

0 

0 . 000 

0.350 

0.350 

0.700 

0.000 

±30 

+0.109 

0.318 

0.318 

0.636 

0.000 

±60 

+0. 182 

0.229 

0.229 

0.458 

0.000 

±90 

+0.200 

0. 102 

0. 102 

0.204 

0.000 

±96* 

+0. 197 

0.075 

0.075 

0.150 

0.000 

30 

-112* 

0. 183 

0.134 

0.112 

0.246 

0.048 

-90 

0.115 

0.182 

0.075 

0.257 

-0.152 

-60 

0.068 

0.240 

0.100 

0.340 

-0.241 

-30 

0.000 

0.318 

0.243 

0.561 

-0.198 

0 

-0.109 

0.318 

0.318 

0.636 

-0.200 

30 

-0.210 

0.248 

0.325 

0.573 

-0.192 

60 

-0.309 

0. 104 

0.251 

0.355 

-0.229 

68* 

-0.440 

0.022 

0.114 

0.136 

-0.387 

60 

-90* 

0.083 

0.269 

0.163 

0.432 

-0.047 

-60 

0.000 

0.249 

0.088 

0.337 

-0.291 

-30 

-0.068 

0.240 

0. 1 00 

0.340 

-0.378 

0 

-0.182 

0.229 

0.229 

0.458 

-0.346 

30 

-0.309 

0.104 

0.251 

0.355 

-0.358 

38* 

-0.449 

0.020 

0.109 

0.129 

-0.483 

90 

-34*/2 14* 

+0.116 

0.196 

0.089 

0.285 

+0.362 

-30/210 

+0.115 

0. 182 

0.075 

0.257 

+0.385 

-10/190 

+0.150 

0.134 

0.075 

0.209 

+0.414 

0/180  * 

+0.200 

0.102 

0.102 

0.204 

+0.400 

6  / 174* 

+0.267 

0.067 

0.123 

0. 190 

+0.411 

TABLE 

IV.  Coupled 

resonator  pe 

rformance  for  case 

of  inhomogeneous 

medium. 

Input  conditions 

given  in  Table  II. 

Asterisk 

denotes 

end  of  stable  operating  regime 

to  within  2° 

(Continued) . 

(a)  $  Fixed 

and  tji 

Varied 

180 
—  4- 

180  J 
~  * 

r 

X1 

X2 

Vl2 

A2 

0 

0 

0.000 

0.350 

0.350 

0.700 

0.000 

±30 

+0. 109 

0.318 

0.318 

1.636 

TO. 200 

±60 

+0. 182 

0.229 

0.229 

0.458 

TO. 346 

±90* 

+0. 200 

0. 102 

0.102 

0.204 

TO. 400 

30 

-92* 

0.122 

0.180 

0.080 

0.260 

0.376 

-90 

0.115 

0.  182 

0.075 

0.257 

0.385 

-60 

0.068 

0.240 

0.100 

0.340 

0.378 

-30 

0.000 

0.318 

0.243 

0.561 

0. 198 

0 

-0.109 

0.318 

0.318 

0.638 

0.000 

30 

-0.210 

0.248 

0.325 

0.573 

-0.192 

60 

-0.309 

0.104 

0.251 

0.355 

-0.358 

68 

-0.440 

0.022 

0.114 

0.136 

-0.494 

60 

-142* 

0.199 

0.137 

0.135 

0.272 

0.122 

-140 

0.197 

0.145 

0.140 

0.285 

0.124 

-130* 

0.182 

0. 183 

0.161 

0.344 

0.141 

-74* 

0.036 

0.269 

0.130 

0.399 

0.242 

-60 

0.000 

0.249 

0.088 

0.337 

0.291 

-30 

-0.068 

0.240 

0. 100 

0.340 

0.241 

0 

-0.182 

0.229 

0.229 

0.458 

0.000 

30 

-0.309 

0. 104 

0.251 

0.355 

-0.229 

38* 

-0.449 

0.020 

0. 109 

0.129 

-0.419 

90 

-60*/240* 

?0. 083 

0.269 

0.163 

0.432 

±0.047 

-50/230 

TO. 103 

0.245 

0.139 

0.384 

±0.070 

-40/220 

TO.  114 

0.216 

0.109 

0.325 

±0. 102 

-30/210 

TO.  115 

0. 182 

0.075 

0.257 

±0.152 

-10/190 

TO. 150 

0.134 

0.075 

0.209 

±0.111 

0/180 

8/172* 

TO. 200 

0. 102 

0.102 

0.204 

0.000 

TO. 340 

0.041 

0.126 

0. 167 

TO. 226 

stable  operating  range  for  ip  .  In  all  cases  the  range  of  xp  for  stable 
operation  is  approximately  ±  90°  which  corresponds  to  maintaining  control  of  h 
to  within  half  a  wavelength. 


B.  FIXED  AND  $  VARIED 

The  present  computations  correspond  to  keeping  L^(Aj  =  0)  and  h(iji) 
fixed  while  varying  <j>)  •  Peak  power  occurs  at  $  =  i|>  =  0  and  for  #  0, 

has  a  local  maximum  at  <J>  =  0.  The  results  for  =  0  agree  with  Eqs.  (23) 
through  (26)  and  indicate  a  stable  range 


_Irl)1/2/Ll  L2 


X  /  2 


<-kI) 


1/2 


( 34 ) 


The  range  of  and  therefore  of  (L^  -  L^)/\  for  stable  operation  decreases 
as  ip  increases.  Tables  Illb  and  IVb  indicate  that  when  ip  t  0  ,  the  solution, 
in  some  cases,  is  a  multivalued  function  of  A^»  The  impact  on  coupled  laser 
performance  has  not  been  pursued. 

C.  GAIN  MEDIUM  EFFECTS 

The  performance  of  coupled  lasers  with  homogeneous  and  inhomogeneous 
gain  mediums  are  compared  in  Tables  III  and  IV,  respectively.  Initial 
conditions  of  Table  II  are  used.  It  is  seen  that  for  these  initial  conditions 
the  range  of  stable  operation  and  the  dependence  of  output  power  on  <}>  and  tj>  is 
insensitive  to  the  nature  of  the  gain  medium.  [The  factor  two  difference  in  I 
between  Tables  III  and  IV  Is  a  gain  medium  property  and  is  not  signi f leant .  ] 
However,  this  good  agreement  may  not  apply  for  initial  conditions  which  yield 


V.  MULTILINE  COUPLED  RESONATOR  PERFORMANCE 


The  previous  results  have  assumed  a  single  lasing  transition.  We  now 

consider  the  effect  of  multiline  operation.  To  achieve  good  beam  quality  it 

is  necessary  to  operate  under  conditions  such  that  <J>  =  0  as  noted  in 

Eq.  (15).  The  stability  of  this  configuration  requires  consideration. 

We  first  consider  a  homogeneous  gain  medium  with  a  single  longitudinal 

mode  for  each  laser  line.  Let  X  and  X  denote  the  smallest  and  largest  wave- 

p  q 

lengths,  in  the  multiline  coupled  resonator,  and  let  X  denote  a  mean  value. 

If  a  lossy  etalon  is  assumed,  the  corresponding  values  of  \Ji  are  if  =  2tt  h/X 

P  P 

and  ili  =  2ir  h/ X  .  Consider  a  typical  device  wherein  (X  -  X  )/X  =  0(0.1)  and 

q  q  q  p 

h  =  0(lm).  It  follows  that 


X( 


(35a) 


=  0(105) 


(35b) 


It  may  be  concluded  that  \p  and  ill  are  uncorrelated  and  the  value  of  ill . 

P  q 


for 


each  wavelength,  is  random  within  the  interval— it  <  iji  <  it.  Since  stable 


operation  is  limited  to  the  interval  -tt/2  <  i|i  <  tt/2  (e.g.,  Fig.  4)  it  follows 
that  for  any  fixed  value  of  h,  and  a  "lossy"  etalon,  only  about  half  of  the 


possible  laser  wavelengths  can  achieve  steady-state  operation. 


In  the  case  of  an  inhomogeneously  broadened  medium,  there  may  be  many 

longitudinal  modes  associated  with  each  laser  line.  Let  X  and  X  denote  the 

P  q 

wavelength  of  adjacent  longitudinal  modes.  Neglect  mode  pulling  and  pushing 

and  assume  (T/R)122  <<  1  in  Eq .  (7b).  It  follows  that  X  -  X  =  X‘/2L  so  that 

q  P 


=  it  h/L 


(36) 


'b  -  'l> 

p  q 

In  the  limit  h/L  +  0,  2,  4,...  all  the  longitudinal  modes,  corresponding  to  a 
given  laser  line,  have  the  same  degree  of  stability  or  instability.  In  the 
limit  h/L  1,  3,  5,...,  in  general,  one  longitudinal  mode  will  be  stable 
while  the  adjacent  mode  will  be  unstable.  In  the  latter  case,  half  of  the 
longitudinal  modes,  for  a  given  line,  are  stable  while  the  remaining  modes  are 
unstable. 

The  net  power  and  stability  of  the  output  from  multiline  coupled  lasers 
will  depend  on  the  competition,  if  any,  between  the  stable  and  unstable 
resonator  modes.  It  is  seen  from  the  coupling  terms  in  Eqs.  (3a)  and  (3b) 
that  when  |$  ±  |  <  tt/2  the  stable  modes  have,  in  effect,  more  gain  than  the 

corresponding  unstable  (i.e.,  uncoupled)  modes.  Hence,  when  both  the  stable 
and  unstable  modes  compete  for  the  same  gain  medium,  the  stable  modes  will 
dominate.  In  the  latter  case,  multiline  coupled  lasers  will  tend  to  provide 
stable  output  on  a  reduced  number  of  lines  with  relatively  little  loss  in  net 
output  power. 
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VI. 


CONCLUDING  REMARKS 


The  present  study  has  considered  coupled  Fabry-Perot  resonators  with 
intercavity  loses  such  that  =  2v  h/\.  In  more  complex  resonator  geometries 
(e.g.,  hole-coupled  and  edge-coupled  unstable  resonators),  the  radiation 
generally  makes  a  single  pass  in  t lie  intercavity  region.  Hence,  a  phase 
change  relation  similar  to  <\>  =  2 tt  h/X,  as  well  as  present  results,  are 
expected  to  be  applicable  for  these  devices. 
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APPENDIX 


ETALON  PERFORMANCE 


Mirrors  and  M2^  may  each  be  etalons.  The  combination  of 

M2p  then,  also  forms  an  etalon.  The  steady-state  performance  of  these 
configurations  is  now  noted. 

I.  SINGLE  ETALON 

A  single  etalon  (e.g.,  mirror  M^)  is  illustrated  in  Fig.  A-la. 
Following  the  development  in  Ref.  3,  the  reflection  and  transmission 
coefficients  can  be  expressed 


where 
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r  =  | ( n' -n) / ( n'+n)  | 

5 1 2  =  4lTh12/Xl2 

and  n,n'  denote  index  of  refraction.  The  coefficient  exp(i5  /2)  in 
Eq.  (A-lb)  is  usually  neglected  but  is  retained  here  to  obtain  phase 
information.  Equations  (A-la)  and  (A-lb)  indicate 


(a)  SINGLE  ETALON  (mirror  M12) 
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12 
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12 


(b)  DOUBLE  ETALON 


Fig.  A-la,  Etalon  configurations 


R12  =  1  "  T12  =  4r  sin  ( 6  j  2 /2 )/ [ ( 1  -  r  T  +  4r"  sin“(6 , 2/2 ) . 


r  ,  .  v  -1  r  l  -  r2  ,612^ 
►  =  (-l)  l - 9  COt(_T“)j 


1  +  r 


t  1 r 1  +  r2  rfi12»1 

►  12  =  tan  L - 2  tan  (“J  J 

l  -  r 


It  follows  that 


*12  5  *U  -  *I2  ‘  *  "/2 


II.  DOUBLE  ETALON 

The  etalon  formed  by  the  combination  of  mirrors  Mj2  and  M2]  is  now 
evaluated.  The  configuration  is  illustrated  in  Fig.  A-lb.  We  assume  that 
scattering  losses  occur  in  the  space  between  the  mirrors.  The  local 
scattering  coefficient  is  taken  to  be  a  constant,  namely 


a  =  -ddnl/dx  =  -2d£nE/dx 
s 


The  reduction  in  electric  field  intensity,  after  one  round  trip,  is 


E  ,,/E  =  e 

P+1  P 


where  a  =  a  h.  The  reflection  and  transmission  coefficients  for  the  double 
s 

etalon,  with  losses,  are  then 
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t  t  T  e'a/2  ei(4,12  +  i/2) 
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where 


6  =  4rrh/X 


5  =  2$j2  +  5 


Equations  (A-5a)  and  (A-5b)  indicate 


R 12  [ (  l  -e  V  +  4e‘°  sin2(«/2)] 


(l  -  R  e  °)2  +  4e  a  R  sin2(6/2) 


( A-6a) 
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2(1  -  R  )  sin(6/2)cos(6/2) 


e  ( 1  -  R12  e  u)(l  -  e  “)  +  2(1  +  R^)  sin"(6/2) 
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laboratory  operations 


The  Aerospace  Corporation  functions  as  an  architect  -engineer  tor 
national  security  projects,  specializing  In  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  1 nves t i gat  Ions  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  Its 
ability  to  stay  current  with  new  developments.  This  expertise  Is  enhanced  bv 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  Che 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state -spec! f ic  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -field -of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  appl Ications ,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communlcat Iona ; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radlometry,  microwave/mll 1 imeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon,  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability,  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy -induced 
env l ronment  s . 

Space  Sciences  Laboratory:  MagneCospher tc ,  auroral  and  cosmic  ray 
physics,  wave -par t  lcle  Interactions,  magnetospher ic  plasma  waves,  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation,  solar  physics,  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere , 
effects  of  e 1 ec t romagnet 1 c  and  particulate  radiations  on  space  systems;  space 
instruroentat ion. 


